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We present the formalism for the decay of dynamically generated scalar mesons with open- or 
hidden-charm and give results for the decay of Z?s (2317) to yDl plus that of a hidden charm scalar 
meson state predicted by the theory around 3700 MeV decaying into yJ/tp- 
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INTRODUCTION 



The discovery of charmed scalar meson resonances at 
Babar, Belle and CLEO 0, H & H has stimulated a 
fruitful line of research, suggesting that their structure is 
much richer than what one might guess assuming the qq 
picture, which has proved quite successful in other areas 
]{| . Some authors have suggested a qqqq structure [f| 0] 
or mixing between the usual qq structure and four quark 
Q. Also there have been suggestions that these states 
might be molecular states of the pseudoscalar mesons 
SMOEO- Similar to this latter works, but with 
subtle differences that we shall discuss later on, are the 
pictures where these states appear as dynamically gen- 
erated in the context of unitarized chiral perturbation 
theory 0, QJ, d E3- Th e works in [3, QJ, El] rely 
upon Lagrangians based on heavy quark symmetry, while 
the one of [l7| starts from an extension of the SU(3) 
chiral Lagrangian to SU(4) which is largely broken due 
to the implicit vector meson exchange characterizin g th e 
Weinberg- Tomozawa term of the chiral Lagrangian llal . 
Because of this, the terms of the Lagrangian used in [l7| 
are suppressed by the ratio of the light to heavy mass 
squared of the vector mesons. SU(A) Lagrangians with 
covariant derivatives adapted to include weak interac- 
tions are used in . The Lagrangian in pj} contains the 
Lagrangian used in [3, [l6| , which is suited for the study 
of open-charm resonances like the D* (2317), but it also 
contains other terms which allow to study the hidden- 
charm states. The comparison of the results using this 
Lagrangian with another one using a chiral symmetry 
breaking extension to SU(N) of the SU(3) results [20| 
allows to have an idea of the uncertainties in the results. 
The stability of the D* (2317) was confirmed, while the 
A (3700) state could sometimes become a cusp instead 
of a bound state, but experimentally it would lead to a 
bump in the mass distributions in any case. Hence one 
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is confident that the new state should also be found. 

The radiative decay of resonances has been usually sug- 
gested as a test for their nature [2l|, • In the present 
case, we are interested in the decay D* (2317) — > D* s y 
and A(3700) — * yJ/ifj- The first of these decays has been 
evaluated assuming varied structures for the resonance, 
within quark models, vector meson dominance , lig ht cone 
QCD sum rules, etc [H H, HE [H, M, HE HI HUE HI 
and more recently from the point of view of the Z?* (2317) 
as a molecular state [33|. The A (3700) has also been pre- 
dicted in [34| as a molecular state assuming a reasonable 
interaction Lagrangian between the D and D states. The 
radiative decay of this predicted state is reported here for 
the first time. 



II. FORMULATION FOR THE RADIATIVE 
DECAY 

In a picture of the scalar mesons as dynamically gener- 
ated one needs to couple the photon to the meson compo- 
nents of the coupled channels. Here the reactions studied 
are 

c: (23i7)^ 7 ^: 

A(3700) -> jj/ip 

and the technical way to evaluate them is considering the 
loop diagrams of Fig. [1] 

These are the same diagrams used in the evaluation of 
the radiative decay of the D* (2317) as a DK molecule 
in 33]. The differences are that here we do not need the 
wave function of the D* (2317), we only need the cou- 
plings of the resonance to the PP channels which are 
obtained in the study of these resonances as dynamically 
generated in [TtJ • Another difference is that we have 
more channels than DK, and this has some numerical 
effects in the results because of cancellations of terms in 
the most important channels D + K°, D°K + . Finally, al- 
though it provides a small contribution in (33[, there is 
a term where the photon couples to the D*. This term 
involves a transition Z?* (2317) — ► D* s which is allowed 
for a virtual D* as it would be the case here. Such terms 
mix the longitudinal part of the vector meson propagator 
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D *(2317) 




theory, as shown in [33| , but they vanish in the radiative 
decay amplitude, as we show in appendix A. 

The procedure to evaluate the radiative decay followed 
here for the dynamically generated scalar resonances has 
been tested with success in the decays <j) — ► /o (980)7 an d 
4> — > a (980)7 [H, [H, H3| with the Jo and a resonances 
dynamically generated from the interaction of the lowest 
order meson- meson chiral Lagrangian (4lj . The present 
reaction is the time reversal reaction, in the charmed sec- 
tor, of the radiative cj> decay into a scalar and a photon. 
The same ideas presented here are used in the study of 
the radiative decay of the / (980) and a (980), as dy- 
namically generated resonances, into jp and joj in [22|. 

The channels to which the D* (2317) and the A(3700) 
resonances have appreciable couplings in are the fol- 
lowing: 

£> fl * (2317) : D+ K° , D°K+ , D+ V , D+ Vc 
A(3700) : D+ D~ , D° D° , D+ D~ 

We shall demonstrate that, using arguments of gauge 
invariance, we can overcome the evaluation of the dia- 
gram c) of Fig [T] and, as a consequence, we must only 
evaluate the diagrams of Fig [5] for the £>* (2317) and of 
Fig O for the X (3700) 





FIG. 1: Diagrams considered in the evaluation of the radia- 
tive decay D* (2317) — > 7-DJ. P and P are the pseudoscalar 
mesons that couple to the D* (2317). 




D *(2317) 



FIG. 2: Diagrams needed in the evaluation of the D* Q (2317) 
radiative decay. 



with the scalar meson. However, in order to prevent the 
appearance of a pole of a scalar in the vector meson prop- 
agator, in a covariant formalism like the one we use, this 
transition amplitude must vanish for P 2 = rrig, with P 
and ms the total momentum and the mass of the scalar 
meson, respectively [35l l36| . This is discussed in a re- 
lated work on the radiative decay of axial vector mesons 
[37] ] , but we address the problem in detail in appendix A. 
There is also a diagram with the photon attached to the 
scalar resonance (for charged states), which is shown to 
vanish in the amplitude in [33j due to the lorentz condi- 
tion of the vector meson. 

The two diagrams discussed above, together with those 
considered in figure [l] provide a set of gauge invariant 
terms, as shown explicitly in (33j. The two terms with 
the photon coupling to the external particles in the loop 
diagram play a role in the gauge invariant test of the 



Let us proceed to the explicit evaluation of the dia- 
grams. The amplitude of the diagram of Fig [5] a) is read- 
ily evaluated as: 



-il = 



d 4 q 



H)5 



D* (2317)^PP 



(q + K) 2 — m\ + ieq 2 — m\ + ie 
i 



-(-i)eQ 1 ^( 1 )(q + q + Kr 



(Q - q) 2 -ml + ie 
X (+*)^^j^e,(D:)(q-Q + qr\ v (1) 

where mi and mi are the masses of the upper and lower 
pseudoscalar mesons in the loop diagram, eQ\ is the 
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X(3700) 



X(3700) 




Ds d; 

FIG. 3: Diagrams needed in the evaluation of the X(3700) 
radiative decay. 



charge (e > 0) of the upper pseudoscalar meson, My 
the mass of the vector meson and Gy the VPP coupling 
(Gy = 55 MeV). The constant <to* (23i7)^D+if° provides 
the coupling of the resonance to the D + K°. In Table 1 
we show the calculated couplings for the poles from (l7| 
located at 2317.25 MeV and (3718.93-i0.06) MeV. For 
this last one we have neglected its imaginary part. 

The coupling of the vector meson to the pseudoscalars 
is obtained from the 57/(4) generalization (see [l7[) of 
the ordinary SU(3) Lagrangian. 



Cppv = -^TrQd^flV'*) (2) 

with g = —M v Gy/f 2 and the £17(4) matrices <j> and V 1 
given by 



/ J 



_ n . rig 

V2 ^ Ve V12 



K- 
D° 



i _2_ i _2k_ 

V2 \/6 VJ3 



K° 
D+ 
Pi 



K+ 
K° 



-2y i tic 
VE s/V2 



D+ 



D° \ 
D 
D 

-37,. 



Vl2 



"7c / 
12 ' 



D 



*o \ 



(3) 



V2 V6 i/T2 



K *- 



if* 



D* 



V6 



VT2 



D* 



D* 



D 



Sfi 



(4) 



In eq. ([T]) we have chosen to substitute f 2 in the de- 
nominator by /tt/d following the prescription in [l7| of 



associating / w (the pion decay constant, / w =93 MeV) to T = -ig D , {2 zn)^pp e Qi^v A= M ^ V e^D^e^j) 

the light pseudoscalars and /d =165 MeV to the charmed '" v2 /tt/d 

pseudoscalars. f d 4 q 1 1 

All the other diagrams are easily obtained by changing x J (2ir) 4 (q + K) 2 — m\ + ie q 2 — m\ + it 

the couplings of the resonance to the channel and the ^ 

constant Ay to account for the different VPP vertices. x -— g ~ ( f ( t (7) 

The different values of Ay are given in Table HU ^ ' 2 

By using the Lorentz condition for the photon and the UpQn integration of the q varia ble one has the expres- 

vector meson, siom 



e,{D* + W = (5) r = T^e,(V)cAl) (8) 

e u (j)K v = 0, (6) 

and Lorentz covariance provides the most general form 
the amplitude of eq. (Q]) is simplified and we obtain: for as: 
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TABLE I: Couplings of the resonances to the pseudoscalars. 



with 



Resonance 




tin nn [MoVl 


D: (2317) 


D + K' 


-7358 


rs+ nU 
K T D 


-7358 


DTv 


5993 




1541 


X(3700) 


D+D- 


7353 


D~D + 


7353 


DJD- 


6740 


D~D7 


6740 



TABLE II: Coefficients Xv for the coupling of the vector me- 
son to the pseudoscalars. 



Vector Meson 


Channel 


Xv 


Dt 


D + K^ 


1 


K+D u 


-1 




-2/V6 


DTvc 


2/V3 


J/i> 


D+D~ 


-2/V3 


D~D+ 


2/V3 


DTD7 


-2/V3 




2/V3 



a 
b 
c 



(Q - qf - m\ 



(q + k) 2 



(13) 
(14) 
(15) 



Upon a transformation q = q' + Q(l — x) — Ky we are 
left with the integral: 



d y (y + Q{i-x)-K V y 



(2tt) 4 
q' + Q(l 



(q' 2 + s + ie) 3 
c) - Ky) 



(16) 



with s = Q 2 x(l — x) + 2Q.K{\ — x)y — m 2 + (m 2 — m 2 )x, 
which shows that the contribution to the d term comes 
from: 



tfY K^Q u {\-x)y 
(2tt) 4 (q' 2 + s + ie) 3 



(17) 



where two powers of q' have disappeared from the integral 
and hence it is convergent. The q' integral is also readily 
done following the Fcynman formalism: 



t m» = ag n» + bQ»Q v + cQ^K" + dK»Q v 

+ eKW (9) 

The Lorentz condition of eqs. ([5]) and ([6]) removes the 
contributions of the 6, c and e terms, such that only the 
a and d terms contribute. In addition, gauge invariance 
(which is guaranteed when all the terms in fig. [1] are 
accounted for) T^K^O, implies 6 = and a+dQ.K = 
so, 



-dQ.K 



(10) 



such that only one term is needed in the evaluation. We 
choose to evaluate the d term because it is finite and only 
comes from the diagrams of figs. [2]and[3] The procedure 
outlined here has been used before in the evaluation of 
the <j> -► -/K°K decay 0, S|. 

The amplitude T is now easily written as: 



d 4 q' 



1 1 



(2tt) 4 (q 12 + s + ie) 3 (2tt) 4 2 s + ie 
and the d coefficient is readily obtained as: 



(18) 



cl = 



-9D* s0 (2317)^PP e Ql^V 



MyGy yg 
UfD 8^ 



dx 

l] JQ 



(1 - x)y 



(19) 



As mentioned above, one can see, following the same 
procedure, that the diagram in fig. [T]c) only contributes 
to the ag^ v term of eq. © and thus we do not need to 
calculate it. 

The finiteness of the results is also noted in [33[ where 
the wave function is governed by a range parameter A, 
and the results remain finite in the limit of A — > oo. 

One must sum coherently the contribution of each term 
of the diagrams in fig. [2] to the d coefficient and then the 
radiative decay width is given finally by 



T = —d(Q.Kg** v — K ti Q")e ll {y)e v {'y) (11) r = _L_i \K\2(K.Q) 2 \d\ 2 (20) 

87T "^^(2317) 

The evaluation of d is straightforward following the 
Feynman formalism. We write: where \K\ is photon three momentum in the rest frame 

of the D* (2317). 

The decay of the X(3700) proceeds identically through 
the same lines using the appropriate couplings and 



i _ r 1 . r , i 



= 2 / dx dy 3(12) 



abc J Q J [ a + (b- a)x + (c - b)yf masses in the diagrams of fig. [3] 
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TABLE III: Results 



Diagram 


d [fm] 


L [KeV] 




0.01284 


2.518 


D+K° 


-0.00529 


0.427 


D+f) 


-0.00197 


0.059 


Dfv c 


0.00007 


0.000 


Total 


0.00565 


0.488 


TABLE IV: Results 


Diagram 


d [fm] 


L [KeV] 


D+D~ 


-0.00314 


3.709 


D-D+ 


-0.00314 


3.709 


D+D- 


-0.00129 


0.622 


D7D+ 


-0.00129 


0.622 


Total 


-0.00886 


29.481 



ensemble of sets for the input parameters within a phys- 
ical allowed range and calculate the radiative decay for 
each set of parameters in the ensemble. The uncertainties 
in the results are then given by the standard deviation 
from the mean value calculated: 



SiIi(r-ro 2 

N - 1 

JV 



N ^ 



(21) 
(22) 



Since the radiative decay of the Z3* (2317) is very small 
and the uncertainties are of the same order of magnitude, 
we will separately calculate the standard deviation above 
and under the mean value. The parameters will be gen- 
erated within the ranges [l7l ]: 



III. RESULTS 

In Table IIHI we show the results for the d coefficient 
from each term in fig. [51 

As we can see, the largest contribution comes from the 
K + D° intermediate state. The D + K° is smaller than the 
K + D° since it involves two heavy pseudoscalar propaga- 
tors instead of two light ones. Next and weaker than the 
others is the contribution of the D^rj channel, and finally 
the D^rjc channel provides a negligible contribution. 

Note that the contribution from the two charge part- 
ners in the isospin 1=0 DK channel is destructive. Had 
the -D* (2317) been an isospin 1=1 resonance, the relative 
couplings to the two channels would have been opposite, 
making thus a constructive interference and we would 
have obtained a width of 4 KeV instead of 0.488 KeV, 
a factor eight times bigger. Furthermore, because of the 
destructive interference, the effect of the Dfij channel, 
which is quite small by itself, becomes relevant. Indeed, 
if we neglect the channels with the D+ meson, the width 
obtained is T =0.872 KeV a factor 1.8 times bigger than 
when one takes them into account. Then, one can see 
that the consideration of all the coupled channels of the 
approach is quite relevant, which introduces one novel 
element with respect to the ordinary molecular picture 
[33| where only the dominant KD channel is taken into 
account. 

The results for the X(3700) radiative decay are shown 
in Table IIVI We see that this radiative decay is consid- 
erably larger than for the Z?* (2317). In this case all the 
terms add constructively. 

Next we perform an analysis of the uncertainties in 
the results. The fact that we have obtained a very small 
width, because of strong cancellations, indicates that it 
should be rather sensitive to uncertainties in the input 
used for the evaluation. 

To evaluate the uncertainties we will follow the same 
procedure used in [ItJ ■ We will take a random generated 



M v = 2060 ± 52 MeV 
f D = 182 ±36 MeV 
U = 100 ± 15 MeV 
Wr,j o (23i7) = 2316 ±39 MeV 
9d :o (23u)^dk = -6420 ± 1790 MeV 
9d :o (2317)^d sV = 5250 ± 1430 MeV 
9D- (23i7)^D 3 r, c = 1450 ± 470 MeV 

wx(37oo) = 3698 ± 35 MeV 
gx( 3 7oo)^D+D- = 8089 ±3125 MeV 
9x(37oo)^d+d- = 5339 ±2100 MeV 

(23) 

When we do the exercise for N—500 randomly gener- 
ated parameter sets, we obtain 

^(2317) = 0.475tS KeV (24) 

and 

rx(37oo) = 18.45 ± 13.00 KeV (25) 

It is instructive to compare our results with those of 
[H[ for the D* (2317). They are rather similar. In [H 
the results vary from 0.47 KeV in some approximations 
to 1.41 KeV in other approximations. Our uncertainties 
stem from different sources, couplings, masses, etc, but 
the range of values obtained is very similar. A compari- 
son of these results with those of different quark models 
is made in [33[ and we address the reader to Table III of 
this reference. The results obtained with the present pic- 
ture are in general smaller than those obtained in quark 
models or other pictures. The destructive interference 
between the two components of the main isospin chan- 
nel is the main reason for it. Precise experiments on 
this rate should help us understand better the nature of 
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this resonance. As for the X(3700), its search as a peak 
in some reactions would be a first step. The search for 
its radiative decay could follow and, given the large rate 
predicted, the investigation of this decay channel does 
not look particularly difficult, specially when the ratio 
of the radiative decay of the Z3* (2317) — * jD* to the 
D* (2317) — > n°D* has already been measured 



10675058) and Scientific Research Foundation of Liaon- 
ing Education Department (No. 20060490). This re- 
search is in part supported by the Japan ( JSPS) - Spanish 
(CSIC) collaboration agreement. This research is part 
of the EU Integrated Infrastructure Initiative Hadron 
Physics Project under contract number RII3-CT-2004- 
506078. 



IV. CONCLUSIONS 

We have presented here the evaluation of the radiative 
decay of the open and hidden-charm scalar mesons which 
are dynamically generated from the interaction of two 
pseudoscalar mesons. The calculations have been done 
for the D* (2317) open-charm scalar state and for the 
predicted hidden-charm state X(3700), not yet observed. 
We found very different results for the two states. While 
the -D* (2317) decay into jD* has a width of around 0.5 
KeV, the X(3700) has a width into jJ/ip of the order 
of 20 KeV, a factor forty times bigger. One of the rea- 
sons, but not the only one, was the large cancellation 
between the two charge partners of the isospin compo- 
nent of the DK 1=0 state. With the obvious similarities, 
we also found subtle differences between the DK molec- 
ular picture for the D* (2317) state and the dynamically 
generated picture. The latter one, including more chan- 
nels than just the DK, showed sensitive effects from the 
D s rj state, particularly because of the large cancellation 
found between the dominant D + K° and K + D° states. 
Yet, within the theoretical uncertainties, the final result 
obtained in the two pictures are rather similar. 

We also presented a different technical way to evalu- 
ate the amplitudes which makes the formalism simpler 
and shows immediately the finiteness of the results using 
arguments of gauge invariance. 

Concerning the X(3700) state and its radiative decay, 
the large width obtained for the decay into jJ/i/j should 
make its observation easy, in principle, and we also re- 
called that the predictions on this state were rather solid 
so that it should be observed as a bound state or a strong 
cusp, in both of which cases the radiative decay could be 
investigated. The observation of this state with its rela- 
tively large radiative width would provide a boost to the 
idea of the low energy scalar mesons with open-charm 
and some particular hidden-charm scalar states as dy- 
namically generated resonances. 



APPENDIX A: PHOTON COUPLING TO 
EXTERNAL LINES 

In addition to the diagrams considered in figure [TJ one 
should also include two extra diagrams, for the case of 
the decay of charged particles, where the photon couples 
to the external lines, see figured! 




Y 

FIG. 4: Diagrams with the photon conected to the external 
lines. 

These two diagrams are explicitly considered in [33| 
and shown to be relevant in the test of gauge invariance. 
However, we show here that they vanish in the amplitude 
for on-shell L>* (2317) and D* . 

To prove this we first evaluate the loop function in 
figure 0J 
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J(P 2 )P^(D* S ) 



d 4 q 



9d*dk 



(2tt) 4 (P + q) 2 -m\ + ie 



9d' 30 dk 



: (P + 2q)%(D s *) (Al) 



Both diagrams in figure [4] imply vector-scalar mixing, 
which appears througth the longitudinal part of the vec- 
tor meson propagator. Indeed, let us consider the dia- 
grams of figure [5] for the vector meson propagator. 
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a) 



d: 



+ 



The iteration of the last diagram of figure [5] and the 
sum of all these terms leads to a geometrical series which 
renormalizes the longitudinal part of the vector propaga- 
tor and leads to: 



b) 




FIG. 5: Diagrams for the renormalization of the vector meson 
propagator. 



We have: 

W^(P) 
p(«)(p) 



A 



e»(D*)e"{D*) 

p2 _ M 2 + ie 



jLV 



a 

piipv 

P 2 M 2 



pv 

P 2 



1 



P 2 _ M 2 



(A2) 



(A3) 



Where in eq. (|A3|) we separated explicitly the propa- 
gator into its transverse (first) and longitudinal (second) 
components. Analogously, figure [SJd provides a contribu- 
tion to the D* propagator given by: 



iV^(P) 



pfi pv 
p2M2 



pp. pv 
PUP 

pn pis 



P 2 
M 2 



J(P 2 ) 2 - 



m D* (2317) 



p2 _ 2 



P 2 M 2 P 2 - m 2 n , ^-^J(P 2 ) 2 



(A6) 



D* (2317) M 2 



Now comes an important renormalization condition 
which is the physical requirement that the longitudinal 
part of the vector meson pro pag ator does not contain a 
pole of the scalar meson J4J, |45| ■ This condition is only 
fulfilled if 



J(P 2 



t D* (2317) 



0. 



(A7) 



Next we evaluate the two terms in the amplitude of 
figure 0J 



x (-i)J(Q 2 )Q»e v (D* s ) (A8) 

This term is zero because of the lorentz condition on 
the vector mesom Q"e„(Dt) = 0. This was already real- 
ized and used in [33j. 

Next we look at the diagram which contributes to the 
amplitude in figure HJd: 



iT^ = -iJ(P 2 )P%(P> s *). 



(A9) 



P 2_ M 2^;S ~ l ) J ( P ) P l lV D; (2317)(P) 

A 



As we can see, independently of the jVV coupling, the 



ic 



l D* (2317)>> 



which 



(-z)J(P 2 )P^^ 



P 2 - M 2 + ie 



(A4) 



where 2?d* (2317) (P) is the propagator of the scalar par- 
ticle. 

One can see that the presence of P 1 P U1 in eq. (|A4|) 
eliminates the contribution of the transverse part of the 
vector meson propagator, hence, only the longitudinal 
part contributes, and we obtain: 



pfj, pv 



p2 _ 



(A5) 



term T^ b ' is proportional to J(P 2 = 

we have shown before to be zero due to the renormal- 
ization condition of the longitudinal part of the vector 
meson propagator. 

Our procedure to evaluate the amplitude, hence, relies 
upon: 

1) The whole set of diagrams is gauge invariant. 

2) The diagrams of figure|3]do not give contribution to 
the amplitude. 

3) Only the set of diagrams of figure [1] give contribu- 
tion. 

4) Using gauge invariance and the procedure followed 
through eqs. (9-11), only the d term has to be evaluated, 
to which only the diagrams of figureQJi andQ}) contribute. 



'_D S * (2317) 
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